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ABSTRACT 

A method for measuring the wettability of rocks by low-field nuclear magnetic 
resonance is described which has application for hydrocarbon reservoir 
engineering or development, or in civil engineering. The method essentially 
determines the water wet pore surface and the oil wet pore surface when the 
sample is saturated with water and oil, by measurements of relaxation times (T1 , 
T2) of the sample placed in a nuclear magnetic resonance device, previously 
brought to various water or oil saturation states, and calculating the wettability 
index by combination of values the relaxation times obtained for said surfaces. 



a) positive primary drainage of an initially 100 % water-saturated sample 
(curve 1 ) ; 

b) positive imbibition (curve 2) ; 

c) negative imbibition (curve 3) ; 

d) negative drainage (curve 4) ; and 

e) positive secondary drainage (curve 5)[.]; and wherein 

n Swir is the irreducible water saturation discussed below in the P ractical 
ImDlementation example. Sor is the residual oil satu ration also discussed below in 
the Practical implementation example, and Ad and Ai are surface are as discussed 
below and are used in the definition of the wettability indices Wl. 
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The curves represent evolution of Pc as a fu nction of Sw with Swir 



being a specific value of SW. 



[nooil roOOSI T here are various types of devices allowing the curves of 
Fig.1 to be drawn. 

fOOQSI— rOOOSl In a device referred to as a "porous plate" device, notably 
described in U.S. Patent 4,506,542, the porous rock sample containing two fluids 
in the continuous phase is placed in an elongate cell ended at its two opposite 
ends by capillary barriers permeable to a first fluid. This first fluid Is injected under 
pressure through the first membrane and the pressure difference between the 
injection pressure and the pressure of the fluid discharged at the other end is 
measured. The pressures of the two fluids and the capillary pressure Pc are 
constant all along the sample, and the saturation is assumed to be uniform. 
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BACKGROUND OF THE INVENTION 
Field of the Invention 



[0001] The present invention relates to a method for measuring the 
wettability of rock samples by measuring the relaxation time by low-field nuclear 
magnetic resonance. 



Description of the Prior Art 
[0002] There are well-known methods for determining the wettability of 

rocks to the -contained w ater and oil that th o y may contain , comprising carrying out 

rock drainage cycles, that is displacement of fluids intended to decrease the water 

saturation, followed by imbibitionL L - thi o t o rm rolating Imbibition reiates . to a 

displacement of the-fluids allowing an increase in the water saturation (Sw) of the 

rock. The capillary pressure Pc at one point Is defined as the difference at 

equilibrium between the pressure P(oil) of the oil and the pressure P(water) of the 

water. This parameter maUco sonse j susefuLoniy if the two fluids are in the 

continuous phase in the porous medium. For a water wet medium, only the 

positive values mako sen6e,areuseful. On the other hand, when the medium has 

a mixed wettability, the fluids can remain in the continuous phase for the positive 

as well as for the negative capillary pressures (Pc). 

[0003] For an application of this type, a complete capillary pressure 
measuring cycle generally comprises (Fig.1) : 
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{Q0061 r00071 It is also known to carry out progressive-speed centrifugation 
by means of centrifugation devices such as those described for example in French 
Patent 2,772,477 corresponding to USH Patent 6.185,985 or French Patent 
2,763,690, or European Patent 603,040 corresponding to U.S. Patent 5,463,894 
or French Patent 2,798,734 filed by the Assignee. 

[00071 roOOBI T he sample initially saturated with a first liquid (brine for 
example) is placed (Fig. 8) in an elongate vessel or cup A containing a second 
fluid of different density (oil for example). Cup A is fastened to the end of a 
rotating arm B and a centrifugal force is applied thereto so as to study the 
displacements of the fluids in the sahiple during at least two distinct steps. During 
a first drainage step, the assembly is then subjected to a centrifugal force applied 
aiong the length of the vessel so as to exert an expulsion force which tends to 
cause part of the first fluid to flow out. The pressure field created by centrifugation 
is expressed as a function of the density p, of the radius R and of the angular 
velocity w, by the relation : 'A w^* . p (Rmax^ - R^), for each fluid. The pressure of 
the two fluids at the sample outlet is the same and is equal to zero at the outlet. At 
the same time, the second fluid flows into the sample. The two fluids move in the 
sample until a position of equilibrium is reached where the force due to the 
capillary pressure in the pores compensates for the centrifugal force exerted. A 
measuring sonde is placed in the cup on the side of the sample. The sonde can 
be of capacitive type for example and comprise a metal rod insulated by a ceramic 
sheath. The capacity between the rod and the conducting fluid (brine), which is 
proportional to the immersed height, is measured. With this measuring system. 



the measuring accuracy is 1.5% of the pore volume. The sonde detects the 
position of the interface between the two liquids in the cup and transmits the 
measuring signals to a measuring signal acquisition and control automaton E 
including hydraulic liquid circulation means and an acquisition device. 

roOOSl [00091 During the re-imblbitlon stage, the velocity Is decreased so as 
to study the return of the initial fluid into the sample. The local saturations 
measured with this type of device are calculated by an inversion program from the 
total amount of water expelled from the sample. 

10009^ roc 1 01 A ccording to another method, referred to as_a "dynamic" 

method, a sample is placed in an elongate cell having water-permeable 
membranes at its two ends. At a first end, oil under pressure is directly injected 
into the enclosure. Water Is also injected, but this injection is carried out through 
the membrane and at a lower pressure. At the opposite end, the oil is directly 
discharged whereas the water flows out through the terminal membrane. 
Adjustment of the oil and water injection rates allows the capillary pressure to be 
the same at the inlet and at the outlet of the enclosure, which leads to a uniform 
saturation that can be deduced from the fluids balance. The capillary pressure is 
obtained for example by measuring the difference between the pressure of the oil 
and of the water at the enclosure outlet. Such a method is notably described by 
Brown H.W. in "Capillary Pressure Investigations" , Petroleum Transaction AIME. 
vol.192, 1951. Examples of Implementation are for example described in patents 
European Patent 729,022 corresponding to U^^U.S. Patent 5,698,772 or EP 
Patent 974,839 corresponding to U.S. Patent 6,229,31 2 field by the Assignee. 
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rooioi room a method referred to as semi-dynamic method is also known, 
wherein a rock sample imbibed with a first fluid is confined in a closed cell, 
another fluid under pressure is injected at a first end of the enclosure and the 
opposite end is swept by a low-pressure fluid circulated by a pump which carries 
the drained fluid outside. The device measures the pressure and the saturation of 
the sample, the amount of fluid discharged and the electric resistivity of the 
sample. This method Is implemented for example in French Patent 2,708,742 
corresponding to UtStU.S. Patent 5,679,885 filed by the Assignee. 

f 00111 r00121 Once drainage and imbibition curves are established, it is 
well-known to calculate the wettability indices Wl from the surface areas Ad and Ai 
marked by the positive and negative capillary pressure curves, as shown in Fig.1, 

by the relation JF/j^sajw, = Log-j-- ■ 

SUMMARY OF THE INVENTION 

ro0121 r0013l T he method for measuring the wettability of a porous rock 
sample in the presence of water and oil according to the invention comprises 
determining the water wet pore surface and the oil wet pore surface when the 
sample is saturated with water and oil, and calculating the wettability index by 
combination of the values obtained for said surfaces. 

[00131 r00141 T h'^ method finds applications notably for analysis of rocks 
taken from an underground formation containing or likely to contain hydrocarbons. 
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r00141 rO0151 Knowledge of various parameters, and notably of the 
wettability of the rocks, is also useful notably for enhanced recovery of a 
formation, by injection of a fluid under pressure, and when the fluid (liquid or gas) 
best suited for effluents displacement is to be determined by means of preliminary 
tests. 

[Q0151 roOISl T he invention also finds applications in civil engineering for 
formation hydrology in order to evaluate the degree of pollution of fonmatlons for 
example, or in the building trade for testing building materials notably in order to 
determine waterproofing treatments for example. 

rOQ16l r00171 Determination of the water wet pore surface and of the oil wet 
pore surface when the sample is saturated with water and oil is obtained for 
example by means of measurements of relaxation times of the sample placed in a 
nuclear magnetic resonance device. 

fOQ171 100181 T he wettability index is for example determined by the 
relation : 

where SMw is the water wet pore surface and SIVlo is the oil wet pore surface 
when the porous medium Is saturated with water and oil. 

jfOQiei r00191 A ccording to an implementation mode, the wettability index is 
determined by the following operations : 

a) measuring the characteristic relaxation times of the water-saturated sample ; 
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b) measuring the characteristic relaxation times of the oil in the sample in the 
presence of water, in a zone close to saturation (Swir) of the sample ; 

c) measuring the characteristic relaxation times of the water in the sample in the 
presence of oil, in a zone close to residua! saturation (Sor) ; 

d) measuring the relaxation times of the sample in a state where its 100 % oil 
saturation point is reached ; and 

e) combining the measurements of the various relaxation times so as to obtain 
the index. 

r00191 r0020] A ccording to a preferred implementation mode, the relaxation 
times of operations a) to c) are determined after subjecting the sample to 
centrifugation. 

rQ0201 r002n A ccording to a preferred implementation mode, the relaxation 
times of operation d) are determined after forced displacement of the fluids in the 
sample placed in a containment cell. 

r002n r00221 A n oil whose intrinsic relaxation time (Tb) is as great as 
possible and as close as possible to that of water, dodecane for example, is 
preferably selected. 

fQ0221 r00231 T he characteristic relaxation times are for example those 
corrrccpondi ng corresponding t o either the saturation curves maxima, or to mean 
values of the curves. 



7 



rQ0231 r00241 Experience shows that the measurement of Inmr obtained by 
means of the method is as sensitive but that-it-requires much less time and that it 
is applicable to a large number of samples. 

BRIEF DESCRIPTION OF THE DIRAWINGS 

r002n [00251 o ther features and advantages of the method and of the 
device according to the invention will be clear from reading the description 
hereafter of an embodiment given by way of non limitative example, with 
reference to the accompanying drawings wherein : 

r00251 r00261 Fig. 1 shows examples of capillary pressure curves necessary 
for determination of the conventional wettability index Iusbm (United States Bureau 
of Mines); 

rQ0261 r00271 Figs. 3a-2A and 2fe2B respectively show the distribution of 
relaxation times T2 for a water wet rock saturated with water and oil, and a 
representation of the phase distribution (matrix in hatched lines, water in light grey 
and oil in daricer grey); 

rOQ271 r00281 Figs. 5a-3A and 3te3B respectively show the distribution of 
relaxation times T2 for a rock of intermediate wettability saturated with water and 
oil (same central part as in H§«Fe-Fis^2). and a representation of the phase 
distribution (matrix in hatched lines, water in light grey and oil in darker grey); 
r00281 r00291 Fig. 4 shows, on capillary pressure curves, the saturation 
values where the NMR relaxation measurements are peri^ormed; 
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r00291 f00301 Fig. 5 shows an example of disWbution of relaxation times T2 
at the saturation values given in Fig. 4 (from the bottom up : 100 % water, Swi, 
Sorand 100 % dodecane); 

r00301 r003n Fig. 6 shows the relation between index Inmr and the standard 
index Iusbm; 

roosn rOQ321 Fig, l diagrammatically shows a containment cell with a 
permanent heating circuit by fluid circulation; and 

r00321 — r00331 Fig. 8 diagrammatically shows the layout of a centrifugation 
device used for implementing the method. 

DETAILED DESCRIPTION OF THE INVENTION 

[0034] It should be noted that the NMR analysis technique essentially 
applies to an object to be tested a first static magnetic polarization field Bo 
intended for aligning the initially randomly oriented hydrogen protons nuclei in the 
direction of the field, then a second impulsive magnetic field oscillating at the 
Larmor frequency, perpendicular to the first one, created by coils excited by a 
control signal to carry out a nuclear magnetic resonance experiment. When this 
impulsive field stops, the return of the nuclei to their initial state or relaxation 
generates electromagnetic signals (echoes) which are detected and analyzed. 
The presence of physical parameters of the object is determined from the 
amplitude characteristics of these signals. 
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r00291 ro0301 Fig. 5 shows an example of distribution of relaxation times T2 

at the saturation values given in Fig. 4 (from the bottom up : 100 % water, Swi, 
Sor and 100 % dodecane); 

I'OOSOT r003n Fjq. 6 shows the relation between index Inmr and the standard 
index Iusbm; 

r00311 r00321 Fig. 7 diagrammatically shows a containment cell with a 
permanent heating circuit by fluid circulation; and 

r00321 1 00331 Fig. 8 diagrammatically shows the layout of a centrifugation 
device used for implementing the method, 

DETAILED DESCRIPTION OF THE INVENTION 

[0034] !t should be noted that the NMR analysis technique essentially 
applies to an object to be tested a first static magnetic polarization field Bo 
intended for aiigning the initially randomly oriented hydrogen protons nuclei in the 
direction of the field, then a second impulsive magnetic field oscillating at the 
Larmor frequency, perpendicular to the first one, created by coils excited by a 
control signal to carry out a nuclear magnetic resonance experiment. When this 
impulsive field stops, the return of the nuclei to their initial state or relaxation 
generates electromagnetic signals (echoes) which are detected and analyzed. 
The presence of physical parameters of the object is determined from the 
amplitude characteristics of these signals. 
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r00341 r00351 A new wettability index is constructed by combining values of 
the water wet pore surface SMw and of the oil wet pore surface SMo when the 
porous medium is saturated with water and oil. The index can be calculated for 
example by the relation as follows : 

or by the relation : 

r00351 r00361 Quantities SMw and SMo are obtained for example by 
measuring the dominant relaxation time in a low-field nuclear magnetic resonance 
experiment. 

ro0361 [00371 In fact, in such an experiment, the fundamental relation 
connecting the longitudinal Ti or transverse T2 relaxation time to surface S and to 
the volume of this pore is as follows : 

= h /?, 5 — (2) 

where Tib.2b is the relaxation time of the fluid outside the porous medium. 
Basically, this relation comes from the fact that the molecules in the pore undergo 
diffusion motions and interact with the porous wall during the measuring time (the 
typical relaxation time is of the order of 100 ms). The surface interactions are 
designated by coefficient pi.2 referred to as surface relaxivity. A gradient term is 
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discarded that is important when measurements are performed at a high magnetic 
field (> 0.1T). Relation 2 is strictly valid for a pore saturated with a single fluid. In 
general, natural porous media have a pore size and therefore ratio SA/ 
distribution, A relaxation times distribution is thus generally observed, but this 
does not modify the method described here. 

f 00371 r00381 W hen two fluids are present in a pore within the porous 
medium, the same physical diffusion mechanism is valid, but the distribution of the 
two fluids in relation to the pore surface is of crucial importance. For example, 
when the medium is water wet, the water is at the surface and interacts therewith, 
whereas the oil is inside the pore and does not interact with the suri'ace. When the 
distribution of the relaxation times is measured for such a system, the result of 
Fig, 3a-2A is observed. The water has shorter relaxation times than in the case of 
a 100 % water saturation, because the volume of water V has decreased (see 
equation 1). The oil is characterized by a relaxation time that is the same as when 
the measurement is peri^ormed outside the porous medium (Tib.2b) because there 
is no interaction with the pore surface. Besides, it is well-known that interactions at 
the water/oil interface do not lead to a relaxation that is different from that of the 
oil outside the porous medium, provided that the oil considered is light (refined oil 
such as dodecane). Diagrammatically. the phase distribution is shown in Fig,_ 
2b2B. When the same porous medium has a different wettability and is not cleariy 
water wet, the distribution of the relaxation times is modified (Fig.— 3a^). 
Essentially, the characteristic time of the oil will be shorter as a result of an 
interaction of the oil at the pore surface because the water is no longer the closer 
to the surface. The possible distribution of the two fluids is shown in Fig, 3fe_3B. 
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Such an effect has been obtained using standard procedures allowing to 
reproduce the surface properties of rocks from a petroleum reservoir : cleaning, 
then aging in the reservoir oil at the temperature of the reservoir at the irreducible 
water saturation, 

r00381 f00391 T he relation (2) is generalized to a two-phase water-oil 
system. For the water, the relationship is: 

^\w,2w ^lBw,2B\v 

and for the oil : 

= + p^^^^^ (4). 

r00391 r00401 T hus quantities SMw and SIVIo can be determined by 
measuring relaxation times Ti or T2 in the porous medium, relaxation times Tib or 
T2B of the fluids outside the porous medium, liquid volumes V© or Vw present in the 
porous medium. The interaction constants p can be determined separately, but it 
will be seen that only the ratio of these constants is necessary. 

f004Ql — r004n Relations 3 and 4 are introduced in the base relation 1 .which 
results in the relationship: 
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where Sw=VwA/p Is the water saturation, 
So=VoA/p is the oil saturation, and 
Vp the pore volume. 

Cp^EjiL is the ratio of the relaxivities for the water wet (w) and oil wet (o) 

Po 

surfaces. To sinnplify the notation, indices 1 and 2 have been suppressed. 

r00 4 n [00421 T he nnost suitable saturations for measuring quantities SIVlw 
and SMo thus remain to be determined. For reasons linked with the calculation of 
the relaxation time distributions, we-the Irreducible water saturation (Swir, Fig. 4^ 
which is a specific value of SW ) is selected to determine the oil wet surface in the 
presence of water (SMw), and the residual oil saturation (Sor, Figur e Fig. 4 ) to 
determine the water wet surface in the presence of oil. Also, an oil whose intrinsic 
relaxation time (Tb) is as great as possible and as close as possible to water is 
selected. In fact, if Tb is too small, the method will be limited to porous media 
whose ratio SN is high (small pore sizes), or to media whose surface relaxivity is 
high. Dodecane for example is a refined oil which is suited to the proposed 
measurement because its Intrinsic relaxation time (Tbo) is 1 s, close to the 
relaxation time of water (Tbw, about 2.7 s). Crude petroleum oils should generally 
not be used because their relaxation times are too small and they also have a 
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relaxation time distribution wliich considerably hinders analysis. The practical 
details of the calculation are given below. 

r00 4 21 100431 From the distribution of the relaxation times at the different 
saturations (Seejhe example of Fig. 5), the relaxation time corresponding to the 
dominant peak (Tw at Sor, To at Swi) can be easily determined, and these values 
are used in relation 5. It is well-known that determination of the dominant 
relaxation time is robust and depends little on the distribution calculation process. 
It can be noted that the dominant peak corresponds to that of the oil at Swi and of 
the water at Sor. The presence of the second fluid at these different saturations is 
important from a physical point of view, but it disturbs calculation very little. It can 
be shown that whatever the relaxation time of the second fluid present in small 
amount (water at Swi, oil at Sor), the relaxation time of the dominant fluid does not 
fluctuate much. This aspect is important for the calculation robustness. Fig. 5 
shows the distribution of the transverse relaxation times T2. The longitudinal 
relaxation time Ti can also be used, but this measurement is less favourable in 
general because Ti > T2 but Tib = T2B for a given fluid. 

f00 4 31 f 00441 In order to determine the surface relaxivity ratio Cp, the 
dominant relaxation times Twioo and T0100 respectively at the two saturations 
Sw=:100 % are used and So=100 %. The formula used is as follows : 
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r00 4 41 r00451 Fig. 5 illustrates the two measurements performed. Thus, in 
total, the relaxation times distributions have to be measured at four saturations : 
100 % brine, 100 % oil and the extreme saturations Swi and Sor. It is possible to 
therefore use the centrifugation technique to reach, starting from Sw=100%. Swi, 
then Sor. A single rotation stage at maximum speed allows these values to be 
rapidly reached. For measurement at So=100%, a succession of miscible 
displacements by injecting an alcohol, then the refined oil are used. 

Implementation 

roO^SI r00461 For implementing the method, an NMR measuring device 
conventionally comprising {Fig.7) magnets 1 with an air gap in which coils 2 are 
arranged is used. The coils are connected to an electronic box which generates 
the excitation signals (signals In the radio frequency range for example) creating 
an oscillating electric field for acquiring the response of sample S to the excitation 
signals. A device of this type Is used for example In French patent application 
2,823,308. 

Practical implementation example 

100 4 61 — r00471 T he wettability index is obtained by carrying out for example 
the following succession of operaUon with a reservoir sample whose wettability is 
intermediate : 

1. The sample is 100% saturated with reservoir brine and placed in the NMR 
measuring device (as shown in Figs.7-8 for example) to measure the distributions 
of relaxation times Ti, T2 and to deduce parameter Tw100 of equation 6 ; 
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2. It is then piaced in the centrifugation device (as schematized in Fig. 8) and 
centrifuged at maximum speed in a vessel filled with oil until the irreducible water 
saturation Swir is reached ; the distributions of relaxation times Ti, T2 are then 
measured in the NMR device so as to deduce To of equation 5 ; 

3. After being placed again in the centrifugation device, the sample is 
centrifuged at maximum speed until the residual oil saturation Sor is reached ; it is 
then transferred into the NMR device again to measure relaxation times Ti, T2 and 
to deduce Tw ; 

4. The sample is thereafter placed in a containment cell such as those 
described in the aforementioned patents and a water and oil miscible solvent, 
then oil (dodecane for example) is injected until the 100% oil saturation is 
reached. This point being reached, relaxation times Ti, T2 are again measured in 
the NMR device and the value of T0100 required for detemrtination of Cp is 
deduced. 

r00471 r00481 T hus all the parameters for calculation of wettability coefficient 
Inmr are l<nown. 

Comparison with index Iusbm 

[00481 r00491 Comparison between wettability index Iusbm and the new 
index Ini^r shows a good correlation between these two quantities {Figure 6) for 
reservoir rocks. The measurement of Iusbm can therefore be advantageously 
replaced by a measurement of Inmr, which is as sensitive but requires much less 
time and is applicable to a large number of samples. 
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r00491 r0050! A n implementation is described example whiere a NIVIR type 
relaxometry method is used to measure the surface and the volume of the pores. 
This method is however not limitative. Any other analysis method can be used, 
notably cryomicroscopy. 
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